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ABSTRACT 



Context. The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) 
effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped 
asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric 
lightcurves. 

Aims. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the 
aim to detect possible deviations from the constant rotation caused by the YORP effect. 

Methods. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and 
used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as 
a linear change in the rotation rate in time da>/dt. Values of dw/df derived from observations were compared with the values predicted 
by theory. 

Results. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model 
failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger da)/dt = 
(1.4 ± 0.6) x 10~ 8 radd~ 2 (3cr error), which corresponds to a decrease in the rotation period by 4.2msyr~ 1 . The photometry of 
Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could 
only constrain maximum values of |dw/df| < 8 x 10~ 9 rad d~ 2 for Cerberus, and |doj/df| < 3 x 10~ 8 rad d~ 2 for Ra-Shalom. 

Key words, minor planets, asteroids - methods: data analysis - techniques: photometric 
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1. Introduction 

The anisotropic reflection of sunlight and thermal emission of an 
asteroid results in a net torque that modifies the asteroid's spin 
state in the long term. This effect, called Yarkovsky-O'Keefe- 
Radzievskii-Paddack (YORP), is important for the spin-state 
evolution of astero ids in the size range from ~ 1 m to ~ 40 km 
dBottke et al. 1 12006). In particular, YORP can reorient the aster- 
oid's spin axis, slow down or spin up their rotation, and trig- 
ger a tumbling state. YORP has been found to be a key el- 
ement in understanding the peculia r rotation pole distributio n 
of large Koronis-family members dVokrou hlickv et al. 2003), 
the distributio n of small members in moderately young aster- 
oid families dVokrouhlickv et al]|2006h . and the distribution of 
the rotation rate of small m ain-belt asteroids (Pravec et al. 2008; 
Polis hook & B rosch 2009). Thanks to its ability to bring aster- 
oid rotation toward the fission limit, it has been suggested that 
it is a universal formation mechanism for small binary aster- 
oids dScheeresll2007t iPravec & Harris! 120071: J Walsh et al.ll2008l) 
and asteroid pairs ( Vokrouhlicky & Nesvornv 2008; Pravec et al. 
120101) . Although the importance of YORP is now widely recog- 
nized, the YORP effec t has been detected on on ly three asteroids 
so far: (1862) Apo llo (Kaasa lainen et alj|2007h . (54509) YORP 
dLowrv et alj|2007t iTavlor et aUfcOOTlT and (1620) Geographos 
dDurech et alj|2008l) . Here we extend this list with a new YORP 
detection for asteroid (3103) Eger. 

The primary importance of YORP detection for specific bod- 
ies consists in the possibility of validating the theoretical mod- 
eling of YORP itself. This is especially important after several 
studies have demonstrated possible difficulties in the theoreti- 
cal prediction of the YORP effe ct magnitude, such as a n on- 
uniform distribution of density ( Sch eeres & Gaskehj 120081) . a 
sensitive dependence on small-sca le irregularities in the shape 
dStatlerl 120091; IBreiter et al] l2009t) and the degree of thermal 
beaming (Rozi tis & Greenl [20121) . As a result, some cases of 
nondetections might be as important as the po sitive detections 
dScheeres & Gaskelll 120081; IBreiter et all 120091) . These are the 
cases where the simple YORP theory predicts a possible detec- 
tion for a present data set, yet an accurate analysis of the obser- 
vations does not reveal any sign of it. 

Motivated by the discussion above, we present new pho- 
tometric observations of near-Earth asteroids (1865) Cerberus, 
(2100) Ra-Shalom, and (3108) Eger. We derived a shape model 
and spin state for these asteroids and detected the acceleration 
of the rotation rate of Eger, which we believe is produced by 
the YORP effectQ For the two other asteroids we studied - 
(1865) Cerberus and (2100) Ra-Shalom - we set only upper lim- 
its on the change in the rotation rate. While in the Ra-Shalom 
case this observational bound is about the upper limit of the ex- 
pected theoretical YORP value for a body of its size and he- 
liocentric orbit, in the case of Cerberus, the constraint is much 
stronger. Reminiscent of the case of (25143) Itokawa, the ex- 
pected Cerberus YORP value is a factor of ~ 3 higher than its 
observed limit. 



2. Lightcurve inversion and search for the 
rotation-period change 

We used the l ightcu rve inversion method of 
Kaasalain en & Torppal (12001 1) and iKaasalainen et al.l d2001l) to 
derive asteroids' shape, sidereal rotation period, and spin axis 
direc tion from the ob s erved lightcurves. Following the method 
of IKaasalaine n et al.l (|2003), we assumed that the rotation 
rate at changes linearly with time as a>(t) = a>o + vt. Both 
wo and u = dw/df were free parameters of the optimization. 
We analyzed the difference (measured by the x 1 values of the 
goodness-of-fit) between a constant-period model (v = 0) and 
a model with v + 0. If this difference is small, the data can be 
fitted with the constant-period model and no deviation from the 
uniform rotation is detected. If, on the other hand, the fit for 
v + is significantly better than for v — 0, the change in the 
rotation rate is detected. 

We found a significant change in the rotation period only for 
asteroid Eger. For Cerberus and Ra-Shalom, the constant-period 
model fitted the available data well, and the free parameter v did 
not improve the fit significantly. To estimate the maximum value 
of |dw/df| that still agrees with observations, we increased or 
decreased u and found the value ofy 2 , for which the discrepancy 
between the data and the model became significant. 

2.1. (3103) Eger 



Kaasalain en et al.l (I2002I) have derived a preliminary shape 
model of Eger from 13 lightcurves observe d between 1 9 86 an d 
1997. These l ightcu rv es were obtain e d by Wisniewskil (H"98 



Wisniewskil 
d 19941) . and 
in the Uppsala 



19911). IVelichko et alJ d 19921) . Ide Sanctis et all 



Pravec et alJ Tl998). The data were archived 
Asteroid Photometric Catalog (UAPC) 



Send offprint requests to: J. Durech 

1 The whole lightcurve dataset, shape model parameters, 
and other details are available from the DAMIT website 
http : / /astro . tro j a . mf f . cuni . cz/pro jects/asteroids3D 
(see also Durech et al. 2010). 



dLagerkvist et al.l l200ll) . We observed Eger during its seven 
apparitions between 1997 and 2012. The aspect data for these 
new photometric observations are listed in Table Q] A detailed 
description of the reduction and m easurement procedures at th e 
Wise Observatory can be found in iPolishook & Broschl d2008l) . 
CCD observations and data processing at Kharkiv, Simeiz, 
Maidanak, Lisnyky, and Kitab observatories wer e done in the 
standa rd way, and the details can be found in iKruglv et al.l 
(2002). It is particularly important that during three apparitions 
- 1997, 2007, and 2009 - our observations were not taken 
at the Earth close encounters of Eger, thus providing data at 
new viewing geometries. This is because Eger is presently 
close to the 3/5 exte rior mean motion resonance with the Earth 
(Mil ani et al.1 [l989), and so the viewing geometry at close 
approaches repeats. 

The photometric data set consists of 70 lightcurves cover- 
ing more than twenty years. Some of them were observed at 
unusually high phase angles larger than 75 deg. The usual ap- 
proach of lightcurve inversion - using the convex shape model, 
together with the combination of Lommel-Seeliger and Lambert 
scattering laws - did not succeed in fitting these high-phase 
lightcurves. The large discrepancy between the data and the 
model was an indication of a significant nonconvexity of the 
shape of Eger dDurech & Kaasalainen! 120031) . When we used 
Hapke's scattering model, the fit improved but was still not 
good. However, we obtained a good fit when modeling the shape 
of Eger as a general nonconvex body using the approach of 
IKaasalainen & , Viikinko"sklld2012l) . 

We derived a unique solution for the sidereal rotation period 
P, the change in the rotation rate v, the pole direction A,/3 in 
ecliptic coordinates, and the shape. The best-fit parameters are 
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Date r A a A j3 Obs. 
[AU] [AU] [deg] [deg] [deg] 



1996 07 20.0 


1.156 


0.194 


40.3 


343.6 


12.7 


Si 


1997 02 05.0 


1.459 


0.483 


9.8 


145.9 


11.0 


Kh 


2001 06 24.9 


1.294 


0.422 


41.3 


326.1 


26.8 


Si 


2002 02 16.8 


1.515 


0.544 


11.8 


139.0 


15.9 


On 


2006 06 28.9 


1.269 


0.388 


42.3 


330.3 


25.0 


Kh 


2006 06 30.0 


1.263 


0.378 


42.3 


330.9 


24.8 


Kh 


2006 06 30.9 


1.258 


0.368 


42.4 


331.5 


24.5 


Kh 


2006 07 26.0 


1.122 


0.163 


46.2 


355.6 


1.1 


Kh 


2007 02 10.3 


1.486 


0.510 


9.4 


144.7 


13.8 


SRO 


2007 02 12.3 


1.496 


0.520 


9.6 


143.4 


14.7 


SRO 


2007 02 17.3 


1.520 


0.550 


11.8 


140.3 


16.5 


GM 


2007 02 19.0 


1.527 


0.560 


12.7 


139.5 


17.0 


Si 


2009 03 22.0 


1.902 


1.113 


24.1 


221.3 


35.4 


Wl 


2009 03 29.0 


1.903 


1.076 


22.6 


219.3 


37.3 


Wl 


2009 04 15.9 


1.898 


1.023 


20.3 


211.4 


40.6 


Si 


2009 05 17.9 


1.867 


1.076 


25.7 


196.9 


38.3 


Wl 


2009 05 18.8 


1.866 


1.080 


25.9 


196.6 


38.1 


Wl 


2009 05 24.8 


1.857 


1.107 


27.5 


195.3 


36.8 


Wl 


2011 06 01.0 


1.410 


0.688 


42.3 


317.0 


26.4 


Ab 


2011 06 04.0 


1.395 


0.654 


42.6 


318.7 


26.4 


On 


2011 06 04.0 


1.394 


0.654 


42.6 


318.7 


26.4 


W2 


2011 06 05.0 


1.390 


0.644 


42.6 


319.2 


26.4 


Li 


2011 06 05.0 


1.389 


0.643 


42.6 


319.2 


26.4 


On 


2011 06 05.9 


1.384 


0.633 


42.7 


319.8 


26.3 


W2 


2011 06 06.0 


1.384 


0.632 


42.7 


319.8 


26.3 


On 


2011 06 09.0 


1.369 


0.600 


42.9 


321.4 


26.2 


Li 


2011 06 10.0 


1.363 


0.589 


43.0 


322.0 


26.2 


Ab 


2011 06 10.9 


1.358 


0.579 


43.1 


322.6 


26.1 


W2 


2011 06 25.0 


1.282 


0.434 


44.1 


331.2 


24.2 


On 


2011 06 27.0 


1.272 


0.414 


44.2 


332.5 


23.7 


On 


2011 06 28.0 


1.266 


0.404 


44.3 


333.3 


23.5 


On 


2011 07 02.3 


1.242 


0.363 


44.7 


336.4 


22.0 


On 


2011 07 08.0 


1.211 


0.311 


45.4 


341.2 


19.3 


W2 


2011 07 12.9 


1.184 


0.269 


46.2 


345.8 


15.8 


Kh 


2011 07 22.9 


1.130 


0.196 


50.4 


358.4 


3.7 


Kt 


2011 07 22.9 


1.130 


0.196 


50.4 


358.5 


3.7 


Ab 


2011 07 23.0 


1.129 


0.196 


50.5 


358.6 


3.5 


On 


2011 07 25.0 


1.119 


0.185 


52.0 


1.8 


0.1 


Ab 


2011 07 27.0 


1.108 


0.175 


54.0 


5.4 


-3.9 


Ab 


2011 07 27.3 


1.106 


0.173 


54.4 


6.1 


-4.7 


On 


2011 08 13.4 


1.023 


0.171 


82.0 


52.4 


-42.8 


On 


2011 08 14.4 


1.018 


0.175 


83.4 


55.6 


-44.1 


On 


2011 10 25.3 


0.953 


0.581 


76.6 


147.1 


-33.5 


PR 


2011 10 26.3 


0.956 


0.584 


76.2 


147.7 


-33.2 


PR 


2011 10 27.3 


0.959 


0.587 


75.8 


148.2 


-32.8 


PR 


2011 10 29.3 


0.966 


0.592 


74.9 


149.3 


-32.2 


PR 


2011 10 30.3 


0.970 


0.594 


74.5 


149.8 


-31.8 


PR 


2011 11 01.3 


0.977 


0.598 


73.7 


150.8 


-31.1 


PR 


2011 11 02.3 


0.981 


0.600 


73.3 


151.3 


-30.8 


PR 


2011 11 03.3 


0.985 


0.602 


72.9 


151.8 


-30.4 


PR 


2011 12 09.1 


1.160 


0.587 


58.1 


164.9 


-17.8 


W2 


2011 12 30.0 


1.274 


0.530 


46.1 


166.4 


-8.7 


Ab 


2012 01 30.3 


1.439 


0.488 


17.6 


154.4 


8.9 


PR 



Table 1. Aspect data for new observations of Eger. The table 
lists the asteroid's distance from the Sun r and from the Earth 
A, the solar phase angle a, the geocentric ecliptic coordinates of 
the asteroid (A,/3), and the observatory (Wl - Wise Observatory, 
1 m; W2 - Wise Observatory, 46 cm; Kh - Kharkiv Observatory, 
70 cm; Si - Simeiz, Crimean Astronomical Observatory, 1 m; 
GM - Goat Mountain Astronomical Research Station, 35 cm; 
SRO - Sonoita Research Observatory, 35 cm; On - Ondfejov 
Observatory, 65 cm; Li - Lisnyky, Kiev University Observatory, 
70 cm; Ab - Abastumani Astrophysical Observatory, 1.25 m; Ki 
- Kitab Observatory, 40 cm; PR - PROMPT, 45 cm). 



Relative x 2 color map for the pole direction of (31 03) Eger 




30 80 130 180 230 280 330 
Pole longitude X 



Fig. 1. Statistical quality of Eger pole solutions given in sinu- 
soidal projection of the sky in ecliptic coordinates. The grade of 
shading and the scale bar on the right indicate the x 2 of the fit 
divided by the number of data points. The formally best-fit case 
with A = 226° and /3 = -70° is marked with a full circle. The 
solid line shows the contour with the^ 2 value 6% higher than the 
minimum value of the best-fit solution. It represents our region 
of admissible solutions. 



as follows: A = 226° + 15°, /3 = -10° + 4°, P = 5.710156 + 
0.000007 h (for JD 2446617.0), and v = (1.4+0.6)xl0- 8 radd- 2 . 
These parameters and their uncertainties were derived using the 
weights of individual lightcurves that corresponded to the in- 
verse of the noise. When all lightcurves have the same weight 
regardless of their quality, the formal best-fit value of u is 
1.2 x 10~ 8 radd~ 2 . The fit to the lightcurves for v = is very 
good with only small phase shifts between the model and the 
data. However, using nonzero values for the v parameter im- 
proves the fit significantly. The uncertainty intervals were esti- 
mated from the increase in x 2 when varying the model parame- 
ters. The reported errors correspond to 3<x for the x 1 distribution 
with ~ 4700 degrees of freedom^ Although the relative uncer- 
tainty of v is large, the zero value is outside the formal uncer- 
tainty interval. For v = 0, the fit is significantly worse (30% 
higher x 1 ) than for the best model. We also note that the an- 
gular difference of ~ 55° between our solution and the best-fit 
solution by Kaasalainen et al. (2002) is significant and arises 
from the limited number of lightcurves used to derive the pre- 
vious model. The pole direction is very well constrained and is 
not sensitive to the value of v, the scattering model, or the shape 
parametrization. The 3<x uncertainty region is shown in Fig.Q] 

As an additional check of consistency, we note that our pole 
solution closely matches a constraint t hat was suggested from 
the 1991 and 1996 radar observations dBenner et al.lll997l) . In 
particular, the viewing geometry at the end of July 1996 was 
nearly equatorial for our pole as suggested by this reference. In 
1991 Aug 10, the line of sight was indeed closer to the pole by 
~ 15°. 

The shape model of Eger is shown in Fig. [2] Although some 
of the lightcurves (those published by iPravec et al.l[l9 98) were 
accurately calibrated to standard R or V magnitudes, they did 
not cover the viewing/illumination geometry sufficiently to use 
them for constraining the model. Therefore, we treated all the 
calibrated lightcurves as relative. 

In general, nonconvex models are much less stable than con- 
vex ones with respect to the errors in the data, the level of reg- 



2 The x 1 distribution with v degrees of freedom has mean v and vari- 
ance 2v. For Eger, the number of data points was ^ 4800, the number 
of model parameters = 100, thus the number of degrees of freedom was 
v =a 4700. 
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Fig. 2. The convex (top) and nonconvex (bottom) shape models of Eger shown from equatorial level (left, center) and pole-on 
(right). 



ularization, the scattering model, etc. However, this is the first 
time the nonconvex shape fits the lightcurves that were observed 
at large phase angles significantly better than the convex one. 
The nonconvex model in Fig. [2] is only one of many similar 
models that we obtained using different regularization, scattering 
models, and resolution. Apparently, the model contains many 
details that may lead to misinterpretation. The main concavity 
is likely to be real because it repeats for models with different 
regularization. The other details differ from model to model and 
are instead artifacts of the modeling process. 

The lightcurve fits for models with v = and v + are 
shown in Fig. [3] The difference between the two models is 
clearly visible only for lightcurves from 1986, 1987, and 1996. 
Because both the number of lightcurves and the number of data 
points in the lightcurves increase towards more recent observa- 
tions, there is almost no phase shift between the models for ob- 
servations between 2006-2012. The sign of the phase shift is 
different for different epochs. 

This is further demonstrated in Fig. |4] where the shift in 
phase between the observed and modeled lighcurves is plot- 
ted. Although the scatter in the phase shifts is of minutes, there 
is a clear quadratic trend for the constant period model (top 
panel), which is what we expect if the rotation rate changes 
linearly in time. The best-fit quadratic function corresponds to 
v = 0.8 x 10~ 8 rad d~ 2 . This is different from the value obtained 
by the lightcurve inversion, because with lightcurve inversion, 
the difference between the observed and modeled brightness is 
minimized, not the phase shift between the observed and mod- 
eled lightcurves. In the bottom panel, the phase differencies are 
plotted for the YORP model with v = 1.2 x 10~ 8 radd~ 2 (equal 
weights for all lightcurves). For this value of u, there is no 
quadratic trend, and the points are randomly distributed along 
the zero value. 



However, the period analysis depends critically on the ob- 
servations from 1986/87 and is very sensitive to possible sys- 
tematic errors ; _AJfhough the data came from three d ifferent 
works (IWisniewskilll987l 1 199 it Ide Sanctis et all 1 1994 . which 
should largely eliminate any possible observational errors, there 
is an insufficient check against potential systematic model er- 
rors. Therefore, our detection of the acceleration of Eger's rota- 
tion due to YORP appears plausible, but it will have to be nailed 
down with more observations in upcoming apparitions (2014, 
2016, 2017, 2019, etc.) so that it becomes a robust detection. 
With more observations, the time line will be enlarged and the 
quadratic trend seen in Fig.|4]can be confirmed. 

2.2. (1865) Cerberus 

Cerberus lightcurves are characterized by unusually large am- 
plitude, up to 2.3 mag, which is the largest lightcurve amplitude 
observed so far for any asteroid. The UAPC contains th e pho- 
tometry of Cerberus fro m 1980 (Harris & Younelll989l) 1989 
dWisniewski et al.1 Il997l) . and 1998 dSarneczkv et all Il999). 
while additional p hotometric observations were published by 
Szaboeta n <i2ooih . Until now, any pole solution had not been 
determined. We observed Cerberus during three apparitions in 
1999, 2008, and 2009. The circumstances of these observations 
are listed in Table [2] 

We derived a physical model of Cerberus with the sidereal 
rotation period P = 6.80328 + 0.00001 h and the pole direction 
A = 298° ± 40°,/? = -72° + 10°. The pole direction is well con- 
strained, and the 3<x uncertainty region is shown in Fig. [7] The 
convex shape model (Fig.[5]l is very elongated with semiaxis ra- 
tios a/c ^ 4.5 and b/c ^ 1.5, suggesting that the real shape 
of Cerberus might be bilobed or even consist of a close syn- 
chronous binary. When using calibrated lightcurves, the model 
is even more elongated and flat with a/c ~ 7 and b/c ~ 2. The 
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1996/7/20.0 1996/7/22.0 
1 .6 r— i ' 1 ' ' ■— i 1 .6 r— , 1 ■ r 




_l I I I I I I I I I I I I L 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Phase of rotation Phase of rotation 



Fig. 3. Examples of Eger's photometric data (points) fitted with synthetic lightcurves based on the convex shape model. The solid 
curve corresponds to the best model with the rotation rate accelerated by v — 1 .2 x 10~ 8 rad d 2 , while the dashed curve corresponds 
to the best constant-period model with v = 0. The viewing and illumination geometry is given by the aspect angle 8, the solar aspect 
angle Go, and the solar phase angle a. 



constant-period model fits well with all available lightcurves. 
Introducing the v + parameter into the modeling did not im- 
prove the fit. We estimated the maximum allowed absolute value 
of the period change to be |u| < 8 x 10~ 9 radd~ 2 . If the change 
in the rotation rate were higher, it would be detectable in the 
lightcurve set as a shift between observations and the model of 
about 10°. 



2.3. (2100) Ra-Shalom 

The UAPC contai ns lightcurv e s of Ra-Shalom from 19 78 to 
1997 obtaine d bv [Harris et all d 19921) [Ostro etalJ (11984 . and 
Prav ec et all (fl998). Kaasalainen et al. (2004) used this dataset, 
together with a few Ondfejov lightcurves from Ra-Shalom's ap- 
parition in 2000, to construct a pole and shape model. Because 
their data set was limited, they flagged their model to be only a 
preliminary attempt. Indeed, it is significantly different from the 
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Fig. 5. The convex shape model of Cerberus shown from equatorial level (left, center) and pole-on (right). 
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Fig. 6. Examples of Cerberus' lightcurves fitted with synthetic ones based on the convex shape model. The viewing and illumination 
geometry is given by the aspect angle 9, the solar aspect angle 8q, and the solar phase angle a. 



new one we present here, because we used additional photom- 
etry from 2003 and 2009 listed in Table [3] The derived shape 
model is shown in Fig. [8] and the fit to selected lightcurves in 
Fig. [9] While our model is an improvement, it is still not per- 
fect. In particular, the pole direction is still not well constrained, 
with the formally best solution at A = 313° and /3 = -45°. The 
statistically admissible solutions cover an irregularly-shaped re- 
gion on the sky as shown in Fig. [10] The contour corresponds 
to solutions with an 8% higher value of;^ 2 than the minimum. 
Because the number of degrees of freedom is v ^ 2300, the 
8% increase corresponds to 3<x in the x 2 distribution (see the 
footnote on page [3}. The statistical weight of prograde-rotating 
spin solutions is negligible. All solutions in our formal 3cr region 
correspond to a retrograde sense of rotation for Ra-Shalom. That 
our zone of a dmissible pole solutio ns is a subset of a similar zone 
presented by Shep ard et alj ([20081) (see Fig. 7 in this reference), 
who carefully analyzed radar and other wavelength observations 



of Ra-Shalom, gives us an increased confidence in our solution. 
The sidereal rotation period is P = 19.8201 + 0.0004 h. No devi- 
ation from uniform rotation has been detected so far. The max- 
imum value of the change in the rotation rate was estimated to 
-4x 10~ 8 < v < 2x 10- 8 radd- 2 . 

3. Discussion 

Two of our pole solutions, namely those of Eger and Cerberus, 
correspond to a near-extreme obliquity value: we have s =* 176° 
for Eger and s 178° for Cerberus. Both values have about a 10° 
uncertainty in realistic terms. While the formal obliquity value of 
the best-fit solution for Ra-Shalom is 144°, Fig. [TUl indicates 
that basically all retrograde values are possible, including the 
value near 180° (note location of the orbital south pole in this 
figure), so at least two, possibly even all three, of our asteroids 
belong to the most populated class of near-Earth objects, those 
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Fig. 8. The convex shape model of Ra-Shalom shown from equatorial level (left, center) and pole-on (right). 
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Fig. 9. Examples of Ra-Shalom's lightcurves fitted with synthetic ones based on the convex shape model. The viewing and illumi- 
nation geometry is given by the aspect angle 6, the solar aspect angle 6q, and the solar phase angle a. 



with an extremely high va lue of obliquity (lLa Spina et alJ [2004: 
iKrvs zczvriska et al. 2007|). 

The retrograde sense of rotation for all our targets is con- 
sistent with an inde pendent estimate of t he Yarkovsky acceler- 
ation in their orbits. ICheslev et ail (120081) attempted to estimate 
Yarkovsky effect in orbits of all near-Earth asteroids by includ- 
ing a formal along-track acceleration in their orbital fit and by 
estimating the related parameter da I At, nam ely a secular chang e 
in the orbital semimajor axis a (see also iNugent et al.l 12012). 
Because the diurnal variant of the Yarkovsky effect is dominant 
in the dynamics of near-Earth asteroids, a positive/negative value 
for da / dt value implies a progr ade/retrograde sens e of rotation of 
the asteroid (see, for example. iBotfke etaL 2006b. In the case of 
all three asteroids discussed in this paper, ICheslev et al] (2008) 
give a negative best- fit value for da/df, thus requiring a retro- 
grade sense of their rotation. This is especiall y important for Ra- 
Shalom, for which the previous solution by Kaasala inen et al.l 



(12004]) implies prograde rotation, while here we find a retrograde 
solution. 



3.1. Theoretical YORP strengths 

We now proceed with a theoretical estimation of the YORP 
strength for our targets. In a sense, this is merely a consistency 
check that the observed secular increase in the rotation rate v 
(in the case of Eger) can be interpreted as a YORP effect de- 
tection. This is because a number of parameters, such as exact 
size, bulk density, and surface conductivity, are only weakly con- 
strained, and even knowledge of the large-scale surface-shape 
featur es may still not be enough for precise YORP computa- 



tions JStat ler 2009; iBreiter et al.ll2009t iRozitis & Greer] [2012; 



iGolubov & Kruglvil2012h . 
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Fig. 4. The shift in phase between the observed and modeled 
lightcurves for the constant-period model (top) and the model 
with YORP (bottom), v = 1.2 x l(T 8 radd~ 2 . The points in the 
top panel are fitted with a quadratic function. 



(3103 ) Eger We used the model of iBreiter & Vokrouhlickvl 
d201 lb to estimate the YORP value of v for our shape of Eger. 
We considered para meters of the Hapke scatterin g model of 
the E-type asteroids (IBreiter & Vokrou hlickv 2011) and varied 
the value of geometric albedo betwee n 0.4 and 0.6. This i s be- 
cause the latest sol utions presented by Trilling et al.l (1201 Oh and 
lHarris et"aTI (1201 ll) give a little larger size of D 1.78 km 
for a lower albedo value of p 0.39, while previous so- 
lutions had smaller sizes of D 1.5 km and larger albedo 
(Ben ner et al.ll 1997b . We also varied the bulk density value be- 
tween 2.5 and 3 g/cm 3 and assumed surface thermal conductivit y 
of 0.01 W/m/K. However, as was shown bv lBreiter et al.l (12010). 
for ID thermal models, the resulting YORP value of v does not 
depend on the value of conductivity. We scaled the lightcurve- 
inversion shape model (Fig. [2) to have the same volume as a 
sphere of 1 .78 km size. A straightforward use of our model gives 
values for v between 4.1 x 10~ 8 radd -2 and 8.3 x 10~ 8 radd~ 2 . 
However, we note that the lightcurve inversion technique does 
not apriori constrain the axes of the shape model to coincide 
with the principal axes of the inertia tensor. Obviously a large 
difference would signal a suspicious model. Assuming a homo- 
geneous density distribution, we find there is about a 5° tilt be- 
tween the z-axes of the two systems, which is too small an an- 
gle for the correction to be determined by photometry analysis. 



Table 2. Aspect data for new observations of Cerberus. The 
table lists the asteroid's distance from the Sun r and from the 
Earth A, the solar phase angle a, the geocentric ecliptic coordi- 
nates of the asteroid {A,/3), and the observatory (On - Ondfejov 
Observatory, 65 cm; Kh - Kharkiv Observatory, 70 cm; Si - 
Simeiz, Crimean Astronomical Observatory, 1 m; Le - Leura 
Observatory, 36 cm; Ma - Maidanak Observatory, 1 m; Wi - 
Wise Observatory, 1 m; HP - Observatoire de Haute-Provence, 
1.2 m; Li - Lick Observatory, 1 m). 



Relative % color map for the pole direction of (1865) Cerberus 
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Fig. 7. Statistical quality of Cerberus pole solutions given in si- 
nusoidal projection of the sky in ecliptic coordinates. The grade 
of shading and the scale bar on the right indicate the x 1 °f the 
fit divided by the number of data points. The formally best-fit 
case with A = 298° and /3 = -72° is marked with a full circle. 
The solid line shows the contour with the^ 2 value 10% higher 
than the minimum value of the best-fit solution. It represents our 
region of admissible solutions. 



Assuming thus that the body rotates about its shortest axis of 
the inertia tensor, and preserving the body's shape, we now ob- 
tain v values between 4.4 x 10~ 8 radd~ 2 and 7.5 x 10 -8 radd~ 2 . 
Compared to the observed value v = (1.4 + 0.6) x 10 -8 radd~ 2 , 
our estimates are a factor of three to six ti mes higher, a situ- 
ation similar to th e case of (54509) YORP (lLowrv et all 120071: 
Taylo r et afll2007l) . In the latter case, the lightcurve and radar 
data were not able to sample signal from about 30% of the sur- 
face, which was one of the sources of the difference between 
the observed and computed v values. In the case of Eger, we 
only have knowledge of the large-scale features of the asteroid 
shape - the lightcurve data do not give any information about the 
small-scale irregularities of the shape and the effects of thermal 
beaming. We believe this is the main source of the difference 
between the observed value for v and the one computed from 
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Table 3. Aspect data for new observations of Ra-Shalom. 
The table lists the asteroid's distance from the Sun r and from 
the Earth A, the solar phase angle a, the geocentric ecliptic 
coordinates of the asteroid (A, /J), and the observatory (Si - 
Simeiz, Crimean Astronomical Observatory, 1 m; Wl - Wise 
Observatory, 1 m; W2 - Wise Observatory, 46 cm; Kh - Kharkiv 
Observatory, 70 cm; On - Ondfejov Observatory, 65 cm). 
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Fig. 10. Statistical quality of Ra-Shalom pole solutions given 
in sinusoidal projection of the sky in ecliptic coordinates. The 
grade of shading and the scale bar on the right indicate the x 1 of 
the fit divided by the number of data points. The formally best-fit 
case with A = 313° and /3 = -45° is marked with a full circle. 
The solid line shows the contour with the x 2 value 8% higher 
than the minimum value of the best-fit solution. It represents our 
region of admissible solutions. The pole of the orbital plane cor- 
responding to 180° obliquity, shown by the triangle, is located 
close to this zone. 



It is interesting to note that the photometric data are consis- 
tent with a highly nonconvex, bilobed model, but they might also 
be fitted w ith a very close binary system . U sing the analysis per- 
formed bv lBellerose & Scheeresl (120081) or lScheeresI (120091) . we 
note there are possible stable equilibria of two ellipsoids (or an 
ellipsoid and a sphere) with reasonable density values between 
1 and 3 g/cm 3 that could match a synchronously rotating system 
with Cerberus' observed period of ~ 6.81 h and small separa- 
tions. 



YORP theory. Indeed, iBenner et al. | (119971) provide a hint from 
the analysis of the radar ranging to Eger that these small-scale 
features might actually be very significant. 

(1865) Cerberus The shape of Cerberus derived above sug- 
gests that even its large-scale structures were not accurately de- 
termined by our convex model. Obviously, we could have tried 
to resolve some nonconvexities in the model, too, but the cur- 
rently available set of the photometric observations is not large 
enough to derive a unique model. We thus consider our convex 
model as the current, but certainly preliminary, state-of-the-art 
representation of its shape. We used t he mean Hapke parame- 
ters o f the S-type asteroids in this case dBreiter & Vokrouh lickv 
|201 lb and scaled the shape model to have the same volum e as 
a sphere with the diameter of 1.6 km (Mai nzer et al.ll20~TTh . We 
also varied the geometric albedo value between 0.1 and 0.2. The 
bulk density was 2.5 g/cm 3 . With these values, we obtained the 
YORP-predicted v ^ 2 x 10 _8 radd~ 2 , which was two to three 
times larger than the conservative bound we obtained from the 
observations in Sect. 12.21 Obviously, this difference may be due 
to a number of unconstrained factors, the inaccurate shape model 
first of all. The key to further analysis of the Cerberus case con- 
sists in obtaining additional observations, which are needed not 
only for possible future detection of the YORP effect but also for 
improving the shape model itself. A good opportunity occurred 
in October 2016, when this asteroid could be observed at large 
phase angles, and from September to November 2017. 



(2100) Ra-Shalom In the case of Ra-Shalom, we scaled our 
shape model to be equivalent to a sphere of ~ 2 3 km size 
dShepard et al.1l20^lTrillinretini20Tol:rHarris et al.ll20Tlb and 
took a geometric albedo of 0.1 (Ra-Shalom is an S-type aster- 
oid). We used a bulk density of 2 g/cm 3 and our best-fit for- 
mal solution for the pole. The formal obliquity value s 144° 
tends to minimize the computed effect because YORP has been 
found to have a node v ^ near e ^ 125° for generic 
shape s (jCapek & Vokrouhlickv 2004; IVokrouhlickv & C apek 
2002). With these parameter values, we obtained v - -6 x 
10~ 8 radd~ 2 , a little larger than its upper bound derived from 
the lightcurve observations (Sect. 12.31 . For this value, we had a 
negative value of v = dw/df, implying that, for a given shape 
model, YORP should decelerate the rotation rate of the body. 
This conclusion would seem consistent with the rather long ro- 
tation period of 19.82 h of this asteroid. However, we should 
point out the large current uncertainty in the obliquity of the pole 
solution (Fig. [T0l). For instance, sampling the uncertainty region 
by pole positions along the meridian of the formally best solution 
(i.e., keeping longitude equal to 313°) and taking latitude values 
from -10° to -75°, we obtain v values from 2 x 10~ 8 radd~ 2 to 
-9 x 10~ 8 radd~ 2 . As expected, for obliquities below ~ 125°, 
we obtain a positive value for v; however, for obliquities greater 
than ~ 125°, which includes the formally best-fit solution, the 
value for v is negative. While positive u values still cannot be 
excluded, the majority of pole solutions in the uncertainty inter- 
val shown in Fig.fTOlcorrespond to negative u values. As a result, 
if carefully observed during the next few years (notably in late 
summer 2013 and 2016), Ra-Shalom might become the first as- 
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teroid for which YORP will be found to decrease the rotation 
rate. 
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